34 International Conference on High Energy Physics, Philadelphia, 2008 

CP Violation in Hadronic Penguins at BaBar 
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We present preliminary measurements of time-dependent CP- violation parameters in the decays B° — > u)K%, B° — » 
r]'K°, B° -* tt°K%, B° -> <j>K° s -K { \ and B° -» K+K~K%, which includes the resonant final states 4>K% and / (980)A"g. 
The data sample corresponds to the full BABAR dataset of 467 X 10 6 BB pairs produced at the PEP-II asymmetric- 
energy e + e - collider at the Stanford Linear Accelerator Center. 

"66 : 

Q 1. INTRODUCTION 

O ; 

Measurements of time-dependent CP asymmetries in B° meson decays through b — * ccs amplitudes have provided 
crucial tests of the mechanism of CP violation in the Standard Model (SM) [l| . These amplitudes contain the leading 
6-quark couplings, given by the Cabibbo-Kobayashi-Maskawa Q (CKM) flavor mixing matrix, for kinematically 
\Q [ allowed transitions. Decays to charmless final states such as 4>K° , ir°K°, r/'K , and ujK° are CKM-suppressed 
b — » qqs (q — u, d, s) processes dominated by a single loop (penguin) amplitude. This amplitude has the same 
r— i weak phase (3 = arg(—V c dV* b /VtdV t * b ) of the CKM mixing matrix as that measured in the b — > ccs transition, but 
is sensitive to the possible presence of new heavy particles in the loop Due to the different non-perturbative 
strong-interaction properties of the various penguin decays, the effect of new physics is expected to be channel 
• dependent. 

tIh ■ The CKM phase is accessible experimentally through interference between the direct decay of the B meson to a 
CP eigenstate and B°B° mixing followed by decay to the same final state. This interference is observable through the 
time evolution of the decay. In the present study, we reconstruct one B° from T(4S) — > B°B°, which decays to the 
CP eigenstate u>K®, r/'K®, rj'K®, Tr°K^, 4>Kgir°, or K + K~K® (Bqp)- From the remaining particles in the event we 



> 



\Q , also reconstruct the decay vertex of the other B meson (-B tag ) and identify its flavor. The difference At = tcp — itag 
of the proper decay times tcp and t tag is obtained from the measured distance between the decay vertices of the Bqp 
and Btag and the boost (J3^y = 0.56) of the T(AS) system. In the tt K® analysis we compute At and its uncertainty 
with a geometric fit to the T(4S) — * B°B° system taking into account the reconstructed K® trajectory, the knowledge 
of the average interaction point (IP) [J], and the average B meson lifetime. The distribution of At is given by 

o 

• • ■ e -|At|/r 

F(At) = 1 =f Aw ± (1 - 2w) [-rjfSf sm(Am d At) - C* cos(Am d Ai)], 

where 77/ is the CP eigenvalue of final state /, the upper (lower) sign denotes a decay accompanied by a B° (B°) 
tag, r is the mean B° lifetime, Arrid is the mixing frequency, w is the mistag rate, and Aw = w(B°) — w(B°) is the 
difference in mistag rates for B° and B° tag-side decays. The tagged flavor and mistag parameters w and Aw are 
determined with a neural network based algorithm 

A nonzero value of the parameter Cf would indicate direct CP violation. In these modes we expect Cf = 
and —rjfSf = sin2/3, assuming penguin dominance of the b — > s transition and neglecting other CKM-suppressed 
amplitudes with a different weak phase. However, these CKM-suppressed amplitudes and the color-suppressed tree 
diagram introduce additional weak phases whose contributions may not be negligible @, 0, B 0] • As a consequence, 
the measured Sf may differ from sin2/3 even within the SM. This deviation ASf = Sf — sin2/? is estimated in several 
theoretical approaches: QCD factorization (QCDF) 

sua, 

QCDF with modeled rescattering [llj, soft collinear 
effective theory 12J, and SU(3) symmetry The estimates are channel dependent. Estimates of AS from 

QCDF are in the ranges (0.0,0.2), (-0.03,0.03), and (0.01,0.12) for uK% tj'K , and 7r°if Respectively 0,0,111; 



SU(3) symmetry provides bounds of (-0.05,0.09) for rj'K and (-0.06,0.12) for ir°K° [IJ. Predictions that use 
isospin symmetry to relate several amplitudes, including the I = I B — » Kit amplitude, give an expected value for 
S^ok near 1.0 instead of sin2/3 [l5| . The modification of the CP asymmetry due to the presence of suppressed tree 
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4>(K + K )K° is at 0(0. 01) lg, [l7(, while at higher K + K masses a larger contribution at 



amplitudes in B 
0(0. 1) is possible [18(. 

In these proceedings, we summarize preliminary measurements of time-dependent CP parameters in the aforemen- 
tioned b — > qqs penguin-dominated B° decays. The u)K°, r/'K , tt°K®, and K + K~K® results are updates of previous 
measurements 



iai2c 

are given in Refs. [23 



22j . while the 4>K®Tr results are first measurements. Detailed descriptions of each analysis 



2J, and [25 1 



2. DETECTOR AND DATASET 

The data used in this analysis were collected with the BaBAR detector at the PEP-II asymmetric-energy e + e~ 
storage ring operating at the Stanford Linear Accelerator Center. We analyze the entire BABAR dataset collected at 
the T(4S) resonance, corresponding to an integrated luminosity of 426 fb _1 and (467 ± 5) x 10 6 BB pairs. 

A detailed description of the BaBAR detector can be found elsewhere [2g]. Charged particle (track) momenta arc 
measured with a 5-layer double-sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH) coaxial with 
a 1.5-T superconducting solenoidal magnet. Neutral cluster (photon) positions and energies are measured with an 
electromagnetic calorimeter, which also provides partial K® reconstruction. Charged hadrons are identified with a 
detector of internally reflected Cherenkov light and specific ionization measurements (dE/dx) in the tracking detectors 
(DCH, SVT). Finally, the instrumented flux return of the magnet allows discrimination of muons from pions and 
additional detection of K® mesons. 



3. ANALYSIS TECHNIQUE 

In the rfK® and K + K~K® analyses we reconstruct the K® in the final states ir + ir~ (K® + ) and 7r°7r° (K° ); 
in the other analyses we use only the tt + it~ final state. Other i3-daughter candidates are reconstructed with the 
following decays: tt° -»• 77; 77 -> 77 (?7 77 ); 77 -> Tr+Tr^Tr (773^); -»• 7? 77 7r+7r- (r^ (ry)7r7r ); r( -> ry 37r 7r+7r- (^ (37r)7r7r ); 
rf — » p°7 (?7p 7 ), where p° — > 7r + 7r~; and u> — » ir + ir~Tr°. The five final states used for B° — ► r]'K® are rj', ^ n K^ + ^_ , 
Vp 1 K° + T T _, v' v (3w)^w K w+7r-^ ^( 77 )^-^°o 7r o, and r/' pi K° 0ji0 . For the B° r)'K° channel we reconstruct the 7/ in two 
modes: rf n[rr)mt and i n( ^ >7T . 

After applying loose selection criteria to reduce the dominant continuum e + e _ — > qq (q = u, d, s, c) background, 
we perform an unbinned maximum likelihood (ML) fit to the data to separate signal from background and obtain 
the CP-violation parameters for each decay channel. As input to the ML fit, we use two kinematic variables, an 
event-shape Fisher discriminant, and, in the (dK®, 4>K®tt , and K + K~K® analyses, resonance masses and decay 
angles. 

In all analyses but tt°K® and Tj'K®, we use, as kinematic variables, the beam-energy-substituted mass 



rriES = J (|s + p • pb) 2 /Eq — pg and the energy difference AE = E* B — \E$, where (E 0l p Q ) and (Eb,Pb) are 
the laboratory four-momenta of the T(45) and the B C p candidate, respectively, and the asterisk denotes the T(4S) 
rest frame. In the tt°K® analysis we use mg, the invariant mass of the reconstructed Bcp, and m m ; ss , the invariant 
mass of the -B tag computed from the known beam energy and the measured Bqp momentum with mass of Bqp 



constrained to the nominal B meson mass 27[. In the r)'K® analysis we use only the AE variable because a mass 
constraint on the B meson during the vertex fit leaves tties and AE completely correlated. 

Further discrimination from continuum background is obtained with the combination of four event-shape variables 
in a Fisher discriminant: the angle with respect to the beam axis of the B momentum, the angle with respect to 
the beam axis of the B thrust axis, and the zeroth and second momentum- weighted angular moments Lq and L2, 
defined as Li — J2jPj x l C0S ^jT > where 6j is the angle with respect to the B thrust axis of daughter particle j, 
Pj is its momentum, and the sum excludes the daughters of the B candidate. In the r)' K® analysis we also use the 
continuous output of the flavor tagging algorithm as input to the discriminant. 
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Table I: Preliminary fit results for signal yields and CP parameters. The first errors are statistical and the second are systematic. 
See Sec. [4] for explanation of results. 



Mode 


Sign&l "Yield 


'It J J 






163 ± 18 


°- 55 ±0.29 ± 0-02 


— 0.52±g;|g ± 0.03 


rj'K 

Wk2. 
V'K° L 


2515 ± 69 

A-iJtJiJ -1— JO 

556 ± 38 


0.57 ± 0.08 ± 0.02 
53 + 08 + 09 
0.82 ±0.19 ±0.02 


-0.08 ±0.06 ±0.02 
—0 11 + 06 + 09 
0.09 ±0.14 ±0.02 


*°K% 


556 ± 32 


0.55 ±0.20 ±0.03 


0.13 ±0.13 ±0.03 


Mode 


Signal Yield 


Peff 


A C p 


K + K~ Kg 
4>K° S 

/o(980)*r° 


lull ± 39 
(see text) 
(see text) 


0.52 ± 0.08 ± 0.03 
0.13 ±0.13 ±0.02 
0.15 ±0.13 ±0.03 


0.05 ± 0.09 ± 0.04 
0.14 ±0.19 ±0.02 
0.01 ±0.26 ±0.07 


4>K* (892)° 
^(1430)° 


58 ±3 
172 ± 24 
535 ± 38 
167 ±21 


n 07+0.03 

(see text) 
(see text) 
(see text) 


(see text) 
0.20 ±0.14 ±0.06 
0.01 ±0.06 ±0.03 
-0.08 ±0.12 ±0.04 



The K + K~K2. analysis is designed to account for variations of CP structure and interference over the Dalitz plot. 
We use an isobar model that includes the K + K~ resonances /o(980), 0(1020), Xo(1550), and XcO to extract /3 e ff 
and Aqp (—Cf) from the amplitude and phase information over the Dalitz plot. In the (pKw analysis we measure 27 
parameters that characterize the interference of S, P, and D Kit partial wave amplitudes. We are able to measure 
the single mixing- induced CP- violation parameter (3 e g, which is accessible only through the 4>K^tt° CP eigenstate 
in which we reconstruct just ~ 60 events, by constraining the other 26 parameters, including Aqp for each partial 
wave, with ~ 800 events from the <pK + TT~ self-tagging final state. 



4. RESULTS 

The preliminary fit results for signal event yields and CP parameters are shown in Table HI We report separate 
results for rj' K® and rj'K® in addition to the combined i]'K° results. The K + K~K® results comes from the high-mass, 
non-resonant region of the Dalitz plot (mpc+K- > 1-1 GeV). The total yield in the low-mass region of the Dalitz 
plot {m K+K - < 1.1 GeV), which are mostly <f>K° s and f (980)K° events, is 421 ± 25. The (j)K°TT° yield is the total 
for all partial waves; each 4>Kj yield is the sum of (f>K®ir° and 4>K + ir~ final states events since both contribute to 
the determination of each direct CP parameter Acp. 

All Sf and (3 e g results are consistent with the value of sin2/3 measured in b — > ccs decays 
world averages are sin2/3 = 0.67 ± 0.02 and (3 = 0.37 ± 0.02. All Cf and Aqp results are consistent with zero direct 
CP- violation. These K + K~K® results favor [3 e g ~ 0.37 and rule out at 4.8<r the solution ? — j3 from the trigonometric 
ambiguity of [3 from the measurement of sin 2/3. All results are statistics limited. The dominant systematic uncertainty 
in the rj' K analysis is related to CP structure in the BB background; the dominant systematic uncertainty in the 
K + K~K® analysis is related to the Dalitz model. 



28, 291. The current 



5. CONCLUSIONS 

We present preliminary updates of our measurements of mixing-induced CP-violation parameters in several 
b — > qqs penguin-dominated B decays and the first measurement in the B — > ipK^ir decay. The (f>K®Tr analysis 
demonstrates a novel technique for extracting CP parameters from interfering amplitudes with relatively few signal 
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events. Significant changes to previous analyses include twice as much data for u)K®, 20% more data for other 
analyses, and improved track reconstruction. 
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